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I. INTRODUCTION 


From July 15" to August 17" of 1997, trial runs using 1:24 scale models of a T- 
ACS Auxiliary Crane Ship, a DDG-963 Class Ship, a Commercial Container Ship and 
Lighter Barges were conducted in the seakeeping basin of the David Taylor Model Basin 
at Carderock Division Naval Surface Warfare Center (CDNSWC), Bethesda, Maryland. 
T-ACS Auxiliary Crane ships are ran by the Military Sealift Command Ready Reserve 
Force. These self-sustaining, rapidly-deployable ships support military sea transportation 
needs and are extremely useful in ports that have limited, damaged or undeveloped port 
facilities. The functions of T-ACS ships are to lift and transfer various loads from either 


themselves or adjacent vessels and piers. 


A total of 251 different model-scale trial runs were performed, each running 


approximately 8 minutes, under a matrix of variable conditions: 


Variable Condition: Description: 
Ship Configurations Config I T-ACS in Center 


Container Ship to Port 


Lighter Barges to Starboard 


Config II T-ACS in Center 


DDG-963 Class Ship to Port 


Lighter Barges to Starboard 


Config II] T-ACS in Center 
Lighter Barges to Port 


Ship Heading 45 degree increments from 0 to 360 degrees relative 


to incoming waves 


Sea State Model-scale sea states 3, 3 + swell, 4 and 4 + swell] 
Boom Slew Angle Angles of 0, 45, 90, 270 and 315 degrees 
Rider Block Location Various positions from 0 to 45.3 degrees 


Rider Block Inhaul Angle —- Various angles from 0 to 16 degrees 


Boom Luff Angle Angles of 25, 29.6, 54.5 and 60 degrees 


A composition of wave height, body motion, velocity and acceleration data for the 
vessels was recorded for each run from an array of sensors with their locations 
dependant upon the ship configurations listed above. Results from the 213 
Configuration I trial runs, the most probable ship configuration, were provided for 
frequency response analysis of the raw data to estimate the Full-Scale Response 
Amplitude Operators (RAQOs) of the T-ACS auxiliary crane ship. An RAO is basically a 
measure or ratio of a vessels response to a regular wave of unit amplitude and thus 


defined accordingly: 


RAO = amplitude _ of _ response 


= . (Zubaly,1996,pp. 322) 
amplitude _ of _the _ wave 
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Tupper expresses the method of studying responses in a seaway in the following 


manner: 


“... the apparently random surface of the sea can be represented by the 
summation of a large number of regular sinusoidal waves, each with its 
own length, height, direction and phase.” 


(Tupper, 1996,pp. 104) 


“'.. the response of the ship in such a sea could be taken as the summation 
of its responses to all the individual wave components. Hence the basic 
building block for the general study of motions in a seaway is the response 
to a regular simusoidal wave.” 


(Tupper, 1996,pp. 104) 


These concepts strongly identify the need for multiple trial runs using a matrix of variable 
conditions. The study of responses such as severe rolling in a beam sea or excessive 
pitching and heaving in a head sea provides extremely valuable insight into a vessels 
limitations. Operating or load-handling contraints can then be established based upon 
existing as well as predicted conditions to prevent material damage, downtime and/or, 


even worse, human injury. 
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II. PROBLEM FORMULATION 


A. THEORY 


ie Random Process 


Determining RAQOs begins with periodic sampling of random processes such as a 
changing waveheight or a vessels changing responses to the changing waveheight over an 
interval of time. In this case, each data channel was sampled at a rate of 32.2 Hz over an 
interval of approximately 8 minutes per data run. This sampling produces the raw data 
from which the entire RAO determination process is based upon. But, before the raw 
data can be used, it must be base-lined to remove any offset or bias errors inherent in the 
sampling system to prevent them from contaminating the process. This is done by 
determining the mean of the data points for a given channel and then subtracting this 
mean from each data point. When this step is completed for each data channel, the 


process is ready to continue. 


Ze Energy Spectra and Power Spectral Densities (PSDs) 


The next phase in determining RAOs is to generate a representative energy 
spectrum curve for each data channel comprised of Power Spectral Densities (PSDs) vs. 
wave frequencies (omegas) in which the area under the curve represents the energy 


associated with the changing data. (Papoulias,2000,pp.109) The PSDs for this curve are 
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derived by expansion of the data points per channel into a Fourier series which 
approximates the shape of the energy spectrum at the various omegas. 
(Papoulias,2000,pp. 121) The mean square of the base-lined, raw data record within a 


narrow band of omega centerd at omega js represented as follows: 


me (@) = S(@)A@ where r =data _mean_over _@_ band 
S(@) = energy _ spectra 
Aq@= omega _band 


(Papoulias,2000,pp. 109) 


Thus, as an integrity check along the way which is utilized later in the PSD and RAO 
comparison porgram, the mean square of the whole, base-lined, raw data record should 


equal the area under the entire spectral curve: 
es (@)= Js (w)d@ (Papoulias,2000,pp. 109) 
0 


o Response Amplitude Operators (RAQOs) 


Finally, the RAOs are determined by the ratio of the response PSDs over the 


reference waveheight PSDs as follows: 


where RAO, = response __RAO 


PSD, = response _ PSD 


PSD, = waveheight __ PSD 


(Papoulias,2000, pp. 124) 


B. COMPARISON TOOLS 


li Spreadsheet Results 


Spreadsheet Results are selected PSDs, RAOs and Omegas for specific 
Configuration I trial runs that were initially processed by Mr. Dan Hayden of CODNSWC 
using a spreadsheet program. Only seven data channels were processed per run. These 
channels were the four pitch and roll channels (channels 3-6) as well as the three relative 
TACS/Lighter positions (channels 28-30). As for the specific data runs, only 19, 27, 57, 


67, 390 and 392 were processed. 


The Spreadsheet Results provided a reference to which the RAO Program results 
would be compared. If the RAO Program results were consistant with each of the seven 
processed channels, it was assumed that the remaining channels would be correct as well. 


In order to ensure that the referenced Spreadsheet Results were correct, three additional 


MATLAB methods of computing PSDs were used to validate both the Spreadsheet 


Results and the RAO Program. 


2 MATLAB Methods 


In addition to the FFT method in the RAO Program, three additional methods of 
computing PSDs were utilized as comparison tools to ensure validity of both the 
Spreadsheet Results and the RAO Program: Welch, Periodogram and Multi-taper 
Methods. The PSDs are computed in units of power per radians while the frequencies in 
Hz are later converted to radians per second. Although there are many similarities, each 


method performs the task quite differently. 


a. Welch: [PSDs,F reqs] = pwelch(x,nwindow,noverlap,nfftfs) 
The selectable options when utilizing the Welch method are as described 


as follows: 


Option: Description 
X Base-lined, Raw, Time-based Data Vector 
nwindow Hamming Window Length for Modified Periodograms 


(Default of [] divides data vector into 8 equal-length 


windows with residual data points discarded) 


noverlap Number of Overlapping, Windowed Data Points (Default 
of [] uses 50% overlap from window to window) 
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nfft Length of the Fourier Transform for each Window 
fs Data Sampling Frequency (32.2 Hz) 


(SPT, welch) 


The Welch method produces PSDs by averaging periodograms of overlapping, 
Hamming-windowed sections of the data vector. Hamming windowing reduces possible 
sidelobes in the spectral estimate in order to reveal the presence of weaker components of 
the signal spectrum that may otherwise get hidden. (Marple,1987,pp. 132) Specified- 
length, discrete Fourier transforms (DFTs) of the overlapping windows are computed as 


follows: 


N-l 
X(k+) = Sox(n+)W,” where Wy =e 2eI™ 


n=0 
N = length(x) 
X = DFT 


(SPT, welch) 


Indexing of (n+1) and (k+1) is used since MATLAB vectors run from 1! to N 
instead of 0 to N-1. (SPT,fft) A modified periodogram for each windowed segment is 


then computed: 


ModifiedPeriodogram = S(e’”) = 





(SPT, periodogram) 


The individual periodograms are averaged to produce just one representative 


periodogram. (SPT,periodogram) Finally, the PSDs are computed: 


PSDs = 5 (en) (SPT,periodogram) 
s 


b. Periodogram: [PSDs,Freqs] = periodogram(x,window,nfft,fs) 
The selectable options when utilizing the Periodogram method are as 


follows: 


Option: Description 
X Base-lined, Raw, Time-based Data Vector 
window Window Coefficients for Modified Periodogram of the 


Input Matrix (Used for “m’x’’n” matricies. Default of [] 


implies single column data vector) 
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nfft Length of the Fourier Transform 
fs Data Sampling Frequency (32.2 Hz) 


(SPT, periodogram) 


The Periodogram method produces PSDs in the same general manner as the 
Welch method previously discussed. A modified periodogram is generated which leads 
directly to computation of the PSDs. The differences between the two are that no 
windowing/segmenting of the data vector occurs and only one periodogram is generated. 
(SPT,periodogram) Thus, the Welch method is essentially a refinement of the 


Periodogram method. 


C. Multi-taper: [PSDs,Freqs] = pmtm(x,nw,nfft,fs) 
The selectable options when utilizing the the third, Multi-taper, method 


are as follows: 


Option: Description 
X Base-lined, Raw, Time-based Data Vector 
nw Determines Number of Discrete Prolate Spheroidal 


Sequences (n=2,,-—1) used as Data Tapers for of 


Estimation of PSDs (Default of [] is nw =4) 


1] 


nfft Length of the Fourier Transform 
fs Data Sampling Frequency (32.2 Hz) 


(SPT,pmtm) 


The Multi-taper method is by far the most complex of the three comparison 
methods for production of PSDs. This method combines linear and nonlinear modified 
periodograms to estimate the PSDs by computing each periodogram using a sequence of 
orthogonal tapers or windows in the frequency domain as specified from the discrete 


prolate spheriodal sequences. (SPT,pmtm) 
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II. PROGRAM DEVELOPMENT 


Le GENERAL DESCRIPTION 


The RAO program contains a sequence of events utilizing standard MATLAB 
functions. The validation methods to be discussed later are dependant upon the 
MATLAB Signals Processing Toolbox from which the commands are comprised of 
standard functions within the source code. Each significant event of the RAO program is 


addressed below: 


1. Initial Set-up 


Before processing a data run, both the raw data disc and the output storage disc 
must be placed in their respective designated drives. Designation of these drives must 
conform to the system being used for processing. In order to start the RAO program, the 
operator enters the desired run number to be processed as a MATLAB command in one 


of the following two formats: 


Command Format: Applicable DTMB Runs: 
rao(##) 19 - 99 
rao(###) 101 — 587 
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Once this has been initiated, the program designates which drive the raw, time-based data 
is to be read from as well as which drive the generated jpg format figures and composite 
matrices for the specific run are to be saved to. A floating-point format with 5 digits is 


placed into effect and the following constants are established: 


Constant: Definition: 
ptavg = 50 Number of Sequential PSDs used for Smoothing 
window = 500 Number of PSDs within a Rational Freq Range 


lambda = 24.175 Actual Model Scale (1:24.175) 


freq = 32.2 Raw Data Sampling Frequency in Hz 


Selection of the “ptavg” constant is the unavoidable compromise between display 
resolution and accuracy which will be addressed later. As for the “window” constant, it 
envelops the first 500 smoothed PSDs which encompass a reasonable range of omegas to 
be addressed under full-scale conditions. The remaining constant, “freq”, is the 


frequency at which the raw, time-based data was sampled. 


$ Initial Set-up------------------------------------- 
function raoirun) 


Caeam( “Hs \' J) 
Char(’D:\’ 


sread 
swrite 


format short e 
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DeEavygG «=aeb0> 
waenag@ew = 500; 
tlamoada = 24.175; 
freq See 


2. Establish Proper String Configuration for the Data Run 


Based upon the run number entered by the operator, proper string configurations 
for the run as well as the folder to which results are to be saved are established. The 
program then differentiates between a two digit and a three digit ran number so that the 
run string evolves as “R###” and the folder as “###”. Broad applicability of designators 


such as these are utilized throughout the RAO program. 


% Establish Proper String Configuration for Data Run-- 
sige Soha EY, 


Sil SGheak oR) e. 
Ce =) iGo eial Ona) 
else 
eu =whehar (4k) .- 
end 
srun = numZser (run 
filename = strcat(sl,srun); 


LP Aue aoe 


folder = streabe(sZ,srun) > 
else 
TOLaAeCr ~=—Srun. 
end 
ZF Determine Number of Data Channels to be Processed 


The run number entered enables identification of the number of data channels to 
be processed due to the particular data channel configuration for the run as follows: 


bS 


DTMB runs: Configuration: Data Channels: 


19 — 503 I 46 
504 - 533 I 38 
534 — 587 Il 24 


Each data channel 1s individually labeled per configuration as listed in Appendixes A, B 


and C. 


$ Determine Number of Data Channels to be Processed--- 
Le wae <= 603 


channels = 46; 
elseir sun == 5353 

channels = 38; 
else 

channels = 24; 


end 


4. Establish File-path and Load the Raw, Time-based Data 


Given the run number and drive designation from which the data 1s to be read, 
“filepath” is established as a string in order to execute the MATLAB command “load” to 
load the raw, time-based data matrix for processing. Once the matrix is loaded, its 
dimensions are identified utilizing the MATLAB command “size.” The number of 
columns representing the individual data channels is designated as the variable 
“channels” and the number of rows representing the number of data points collected is 


designated as the variable “n’, both of which will be used throughout the program. 
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% Establish Filepath and Load the Raw, Time-based Data 


filepath = [sread, filename]; 
load (filepath) 

szdata = size(data) ; 
channels = szdata(2); 

n = szdata(1); 


5: Compute both Model-Scale and Full-Scale Omega Ranges 


Before going any further, the Model-Scale omega range is computed based upon 
the sampling frequency, total number of data points and the “window” length designated 
earlier. Dividing the Model-Scale range by the square root of the model scale “lambda” 
produces the Full-Scale omega range. This is supported by the basic assumption for 
presentation of motion data that “Natural periods of motion vary as the square root of the 
linear dimension.” (Tupper,1996,pp.106) In addition, Tupper identifies that “In watching 
model experiments the motion always seems rather ‘rapid’ because of the way period 
changes . Thus, a 1/25 model will pitch and heave in a period only a fifth of the full- 


scale ship. 


@ Compute both Model-Scale and Full-Scale Omega Ranges 
momega = 2*pi*fregq/n* (0:window-1); 
fomega = momega/sqrt (lambda) ; 


6. Full-Scale Scalers, PSD Units and RAO Units Library 


There are four individual sets of Full-Scale scalers, PSD units and RAO units in 
the library associated with the type of data collected by the a particular channel. The 


combinations are as follows: 
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Data Type: Full-Scale Scaler: PSD Units: | RAO Units: 


Wave Height l in dimensionless 
Distance 12/lambda deg deg/in 
Velocity 12/lambda‘(3/2) deg/s deg/in s 
Acceleration 12/lambda‘“(2) deg/s*(2) deg/in s‘(2) 


$ Full-Scale Scalers, PSD Units and RAO Units Library- 


foi = is 

pul = in) 

rul = ‘(dimensionless)’; 
fs2 = 12/lambda; 

pu2 = ‘(deg)’; 

rig = ““(Geg/im) “; 

fs3 = A27Jambaa-(372)° 
pus. =“ Wdec7s) 

rus = (deo 277s)’; 
fs4 = i2/7tambda 2; 

pue = *(deog/s-2)’; 

rudd = *(deg/in s*2)’; 


fe Looping/Processing of Data Channels 


It’s now time to begin looping to process the “channels” sequentially. From the 
first data channel to the Jast, each is assigned a Full-Scale Scaler, “fscaler’; a PSD Units, 
“psdunits”; an RAO Units, “raounits”; and a Channel Title, “chtitle”. These assignments 


will be utilized later for labeling the generated jpg figures. 
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$ Looping/Processing of Data Channels----------------- 
EOr- Gi =. channels 


af Choanneis==46 


ae —— 
fscaler = eorae 
psdunits = pul; 
Paounits = fy. 
cherere = ‘Wave Ht Bow’; 
elseif ch== 
fscaler Sele 
DsGuniCs = puis, 
TAQURT ES sua nate 
chtitle =  Sejaulp<  Sleuiie 


el seam eCh==—26) 


fscaler fs4; 

psdunits = pué4; 

FaQuni Es (= rus: 

chtitle = 'Lgohtr1-PBow TvAcc'; 
else 


end 


elseif senannes s==38 


Ie (es gl 
fscaler Se esis 
psdunits = pul; 
racounles™ = sue 
Gheiele = 'Wave Ht Bow'; 
elseif ch==2 
fscaler SS ec: 
PSaGuni ES seem, 
recuntes m= rua, 
chtitle = SOnisemSsonic’: 


elseif ch==38 
fscaler = sl 
psdunits = pu4; 
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Taouni esa] shu4 
chtitle ‘Lohtr1-PBowTvAcc’; 
else 
end 


elseif channeis==24 


ie Co==) 
fscaler =e 
psdunits = pul; 
raocuni ts meee 
chtitle = ‘Wave Ht Bow’; 
elseif ch==2 
fscaler ease 
DSOQUna eS mou 
Baound nes a= paul - 
chtitle = WSonix Sonre = 


elseif ch==24 


fscaler = fs4; 
Dsemiies. = oud 
raounits = rudd; 
chtitle =. 'Boomrip-Vere Ace > 
else 
end 
end 
8. Base-lining the Raw, Time-based Data 


The actual processing of data begins here by base-lining or removing the mean 
from the raw, time-based data. Otherwise, “Failure to remove large sample means or 
other trends in the data may result in distorted or biased spectral estimates.” 


(Marple, 1987,pp. 132) 
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$ Base-lining the Raw, Time-based Data---------------- 
D4 = SGQatat:) Ga 
>. = x-mean (x); 


2: FFT Computation 


The MATLAB command “fft” returns the Fourier transform of the base-lined, 
column data matrix. If the matrix length is a power of two, the “‘fft” command employs a 
high-speed fast Fourier transform algorithm. But, this is not always the case since the 
number of data points per channel varies from run to run In this situation, an alternate 
mixed-radix algorithm finds the prime factors of the column matrix length then computes 
the discrete Fourier transforms of the shorter sequences. [Either way, the Fourier 


transform of the base-lined data is produced. 


$ FFT Computation 
Y = ine (xX, n) , 


10. Model-Scale PSD Computation 


The Model-Scale PSD is computed by multiplying the elements of the Fourier 


Transform by their complex conjugate and dividing by the total number of elements. 


$ Model-Scale PSD Computation------------------------- 
PSD =i s conga ae 
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11. Moving-window Averaging of Model-Scale PSDs 


This is where the RAO program method differs significantly from the tools in 
MATLAB previously identified by offering a much finer control over the resolution vs. 
display accuracy issue. An envelope initially starts with half the “ptavg’’ window plus 
one PSD points, averages them together and assigns that average as the first, smoothed 
PSD point. The envelope then accepts the next sequential PSD point, averages the now 
larger window and assigns that average as the second, smoothed PSD point. The process 
continues until the envelope reaches it’s maximum size of “ptavg” plus one where the 
latest average is assigned as the smoothed PSD point for the middle of the envelope. 
From here on out, the envelope accepts the next sequential PSD point while dropping the 
first one, thus maintaining a fixed-length, moving window. Refer to the following 


example: 
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Example: Let “ptavg” =4 


[#1 #2 #3 #4 #5 #6 #7 #8 #9 #10] 


[ 1" avg PSD] Window starts at half “ptavg” plus one (i.e. 2 + 1) 
i We ave PSD Window increments in size by one 
[ ci avg PSD ] Max capacity of ptavg + I (i.e. 4+ 1) 
is 4" avg PSD | Window begins shift to the right 
if 5" avg PSD ] Window continues to shift 
[ 6" avg PSD ] Window still shifting 
| aoe eee i 


At this point in the RAO program, the smoothed PSDs are designated as a column 
of data in a composite PSD matrix with the column designation correlating directly to the 
data channel being processed. Thus, the first channel processed becomes the first column 
of PSD data in the composite PSD matrix. As for the chosen window length of 50 plus 


one, tests of trial runs using incrementally larger point average envelopes yielded a 
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compromise of “ptavg” = 50 which offered the best display resolution while maintaining 


accuracy in comparison to the spreadsheet results provided (see figs. 1-3) 


% Moving-window Averaging of Model-Scale PSDs for 
Smoothing 
halfptavg = ptavg/2; 
sum = 0; 
for k = 13peava 
Sums j= sum. + PSD(k, 1)> 
if k > halfptavg 
psd(k-halfptavg,ch) = sum/k; 


else 
end 
end 
pts = window + halfptavg; 
ptavgplusone = ptavg + 1; 
J = i; 
for k = ptavgplusone:pts 
sum = sum + PSD(ptavgt+j,1) - PSD(j,1); 
psd(k-halfptavg,ch) = sum/ptavg; 
5) alae ake 
end 
meansgr = mean(x.*2); 
psdarea = trap(momega,psd(:,ch)); 
psdascaler = meansqr/psdarea; 
psd(-,cH) = psd(:,ch)*psdscaler 
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Figure 3. 


12. Computation of Model-Scale RAOs 


As previously addressed in the theory, RAOs are obtained by taking the 
individual square roots of the desired PSDs divided by the seaway or wave height PSDs. 
With the first data channel for each configuration allocated as the reference wave height, 


the first channel of RAQOs processed will naturally be at unity, a constant dimensionless 


value of one. 


%® Computation of Model-Scale RAOs 


xX = BSG nuChUMy psa) 
mrao(:,ch) Seqrc (xX); 
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13. Computation of Full-Scale RAOs 


Computation of the Full-Scale RAQOs is simply a process of multiplying the 
Model-Scale RAOs by the Full-Scale Scaler. This previously selected scaler is based 


upon the nature of the data sampled. 


weembputalion Of Pull=Scale RAQs—————_—_———— 
fraog(:, ch} = mrao(:,ch) *fiscaler; 


14. Plot Model-Scale PSDs and Save Figure in jpg Format 


A plot of the smoothed, Model-Scale PSDs vs. their associated Model-Scale 
omegas is generated (see fig. 4) then saved as a figure in jpg format in a folder dedicated 
to the particular run and titled accordingly with the run number, channel number and a 


‘“‘p” for Model-Scale PSDs label. 


@ Plot Model-Scale PSDs and Save Figure in jpg Format- 


CLE 

plot (momega (i,:),psd(:,ch)) 
Cpe sel 

SG. = 7, HMUMNZSEr (Gh): 


title([’Model-Scale PSDs for DTMB Run ‘’srun,’ / 
Channel “7 sen, 7: cnhtictie)) 

xlabel(’Omega (rad/s) ’) 

ylabel([’PSD ‘psdunits]) 

saveas(gcf, [swrite, folder, ’\’, filename, ’C’,sch, ’p’], 
wane 
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Figure 4. 


15. Plot Model-Scale RAOs and Save Figure in jpg Format 


A plot of the Model-Scale RAOs vs. their associated Model-Scale omegas is 


generated (see fig. 5) then saved as a figure in jpg format in a folder dedicated to the 


particular DTMB run and titled accordingly with the run number, channel number and an 


‘“m’” for Model-Scale RAQs label. 


$ Plot Model-Scale RAOs and Save Figure in jpg Format- 


(eo) 


ep ha 
pioOE(momeda(il, +), mraodts, ch) ) 


grid 
title([’Model-Scale RAOs for DTMB Run ‘srun,’ / 


Chayne!l *,5¢h;,*s “chtiele]) 
xlabel(’Omega (rad/s)’) 
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ylabel([’RAO ‘raounits]) 
Saveas(gCcl, [swrrte, folder,” \~, filename, ’C’,sch,‘m’ ], 


‘jpg’); 
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Figure 5. 


16. Plot Full-Scale RAQOs and Save Figure in jpg Format 


A plot of the Full-Scale RAQOs vs. their associated Full-Scale omegas is generated 
(see fig. 6) then saved as a figure in jpg format in a folder dedicated to the particular run 


and titled accordingly with the run number, channel number and an “f” for Full-Scale 


RAOs label. 
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$ Plot Full-Scale RAOs and Save Figure in jpg Format— 


Cults 
plot (fomega(1,:),frao(:,ch)) 


Gmia 
title([{'Full-Scale RAOs for DTMB Run '‘'srun,' / 


ehannel ',sch,’s "ehtitle]) 
xlabel('Omega (rad/s)') 


ylabel(['RAO ‘raounits]) 
saveas (gcf, [swrite, folder, '\',filename, 'C',sch, 'f'], 


nog); 


Processing of the first data channel is now complete. The RAO program then 
loops back to event 7 listed above, Looping/Processing of Data Channels, until all data 
channels for the specific run, as determined by the data configuration, have been 


processed to generate the Model-Scale PSDs, Model-Scale RAQOs and Full-Scale RAOs. 
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Figure 6. 
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17. Save Omega Matrices as Last Column of Composite PSD and RAO 


Matrices 


For the purpose of convenient analysis, the single-column omega matrices are 
tacked on as the final columns of their respective PSD or RAO composite matrices. This 
basically eliminates the necessity of the jpg formatted figures by allowing direct 
comparison between any two columns of a composite matrix simply by using a plotting 


function. Thus, PSD vs. PSD, PSD vs. omega, RAO vs. RAO and RAO vs. omega plots 


can readily be produced 


& Save Omega Matrices as Last Column of Composite PSD 
and RAO Matrices 


eh ee oy ce ae 

OSG) Cn) = momega(l1,:)’; 
mrao(:,ch) = momega(l1,:)’; 
imiaqcle ees au = fomega(i1,:)’; 


18. Save Composite PSD and RAO Matrices 


This final event in the sequence saves the composite PSD and RAO matrices 
along with the jpg formatted figures previously saved in the folder dedicated to the 


particular DIMB run. It’s now time to process another DTMB run. 


% Save Composite PSD and RAO Matrices----------------- 

save([swrite, folder, ’\’, filename, ’psd.txt’],’psd’, 
‘-ascii’, ’-tabs’); 

save ([swrite, folder, ‘\’, filename, ’mrao.txt’],’mrao’, 
‘-ascii’, ’-tabs’); 

save ([swrite, folder, ’\’,filename, ’frao.txt’],’frao’, 
‘-ascii’,'’-tabs’); 


Sil 


B. VALIDATION 


Validation of the RAO Program was performed by utilizing a MATLAB 


comparison program, “psdandraocomparisons.m”, written to compare the various 


methods of resultant PSDs; Welch, Periodogram, Multi-taper, Spreadsheet and RAO 


Program. The comparison program is comprised of a sequence of events with each 


explained individually as applied to run 67: 


1. Load Applicable Matrices 


Once the comparison program has been initiated, the following matrices are 


loaded: 


Matrix: 


RO67 


Run6/7/.txt 


RO67psd.txt 


RO67mrao.txt 


RO67frao.txt 


Definition: 

Raw, Time-based Data (All 46 Channels) 
Spreadsheet Model-Scale Omegas, PSDs and RAOs 
RAO Program elel ole PSDs and Omegas 
RAO Program Model-Scale RAOs and Omegas 


RAO Program Full-Scale RAOs and Omegas 


a2 


$ Load Applicable Matrices---------------------------- 
hoad  h06 7; 

HOaG RUNG 7. GC; 

load R067psd. txt; 

fHoO2a AUG /mrao. Exe, 

load ROG7irao. exe; 


2 Base-line the Raw Time-Based Data 


Prior to initiating the MATLAB methods of Welch, Periodogram and Multi-taper, 
the raw, time-based data for both Channel 1: Wave Ht and Channel 3: TACS Roll are 


base-lined for continuity just as it is done in the RAO Program. 


% Base-line the Raw Time-based Data from R067--------- 
x = data(:,1); ® Wave Ht Data 

x = xX - mean(x); 

Z Oata( =, 3) 7 $ TACS Roll Data 

Z = 2% - mean(z); 


3. Compute Wave Ht PSDs using MATLAB Methods 


The Welch, Periodogram and Multi-taper methods are used to compute their 
Wave Ht PSDs and associated Frequencies. The Frequencies will be converted to 
Omegas later. For efficiency of computation, a maximum of 8096 FFTs is designated. 
This happens to be the largest power of two that is less that the total number of data 
points per channel which range from approximately 14,000 to 16,000. Otherwise, 
MATLAB uses a slower approach of seeking out the primes of the total number of data 


points and computing the FFTs according to these segments. 
a2 


%$ Computation of Wave Ht PSDs and Associated Fregs---- 


f[hwelch, fwelich] = pwelch(x, [],[],8096,32.2); 
fhperio, fperio] = Joeriodogram(x, {i xeege 32.2) 
lemme i, fmulei] =Somem(x, |/,e09ons2Z 
meansqrx = mean(x.*2); 

welchareah = trap(wwelch, hwelch) ; 
perioareah = trap(wperio, hperio) ; 
multiareah = Strap (WMmulot. mes.) 
hwelchscaled = hwelch*meansgqrx/welchareah; 
hperioscaled = hperio*meansqrx/perioareah; 
hmultiscaled = hmulti*meansqrx/multiareah; 


4. Compute TACS Roll PSDs using MATLAB Methods 


The same process described above is now performed for the TACS Roll PSDs. 


%$ Computation of TACS Roll PSDs and Associated Freqs-- 


[rwelch, fwelch] = pwetleniz, [/,[)] 76026, 32-2) : 
[rperio, fperio] = periodogram(z, [],8096,32.2); 
femod ei fm ei) = ~onjien(2, [lere0le 5202) 
meansgrzZ = mean(z.*2); 

welcharear = trap(wwelch, rwelch) ; 
perlioarear = trap(wper1o,rperio); 
multiarear =) Crap (wmiuler,rmu es) © 
rwelchscaled = rwelch*meansgqrz/welcharear; 
rperioscaled = rperio*meansgqrz/perioarear; 
rmultiscaled = rmulti*meansgqrz/multiarear; 


5, Convert MATLAB Method Frequencies to Omegas 


Conversion of frequency to omega is done simply by multiplying the frequency 


matrices by 27. 
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$ Convert Frequency to Omega-------------------------- 


wwelch = 2*pi*fwelch; 
wperio =m 72 or espero; 
wmulti = 2 Soskigeee sly. 


6. Compute Areas under the PSD Curves 


As an integrity check along the way, the mean squares of the base-lined data are 
compared to the areas under the PSD curves. At this point, both the Spreadsheet and 
RAO Program results are integrated into the comparison program to ensure continuity 


between the various methods with comparative results as follows: 


a. Wave Ht Comparison 

Method: Result: 
Mean Square of Base-lined Wave Ht Data 0.0124 
Welch Area (rad/s) 0.0124 
Periodogram Area (rad/s) 0.0124 
Multi-taper Area (rad/s) 0.0124 
Spreadsheet Area (rad/s) 0.0126 
RAO Program Area (rad/s) 0.0124 


3D 


From the results listed above for run 67, the Welch, Periodogram, Multi-taper and 
RAO Program methods agree to the nearest 1/1000" when rounded to 5 decimal 


positions. In addition, each is within 1.6% of the Spreadsheet results. 


b. TACS Roll Comparison 


Method: Result: 
Mean Square of Base-lined TACS Roll Data 0.0061 
Welch Area (deg/s) 0.0061 
Periodogram Area (deg/s) 0.0061 
Multi-taper Area (deg/s) 0.0061 
Spreadsheet Area (deg/s) 0.0059 
RAO Program Area (deg/s) 0.0061 


From the results listed above for run 67, the Welch, Periodogram, Multi-taper and 
RAO Program methods agree to the nearest 1/1000" when rounded to 5 decimal 


positions. In addition, each is within 3.3% of the Spreadsheet results. 
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% Compute Areas under the Wave Ht PSD Curves---------- 


welchareah2 = trap(wwelch, hwelchscaled) ; 
perioareah2 = trap(wperio, hperioscaled) ; 
multiareah2 = trap(wmulti, hmultiscaled) ; 
spreadareah = Sram (muno 7 (72) un ow Vege 
raoprogareah = ‘rap ne /psal(274/7) Ree lesaie |) )- 


$ Compare Areas under the Wave Ht PSD Curves against 
meansgrs 

wavehtpsdareacomparison = f[meansgrx welchareah2 
perioareah2 multiareah2 spreadareah raoprogareah] 


% Compute Areas under the TACS Roll PSD Curves-------- 


welcharear2 = trap (wwelch, rwelchscaled) ; 
perioarear2 = trap (wperio, rperioscaled) ; 
multiarear2 = trap(wmulti, rmultiscaled) ; 
spreadarear =""Grap(rune/ (2,2) 2unoy | ae, 
raoprogarear = trap(R067psd(:,47),R067psd(:,3)); 


$ Compare Areas under the TACS Roll PSD Curves against 
meqnsgrs 

rollpsdareacomparison = [meansqrz welcharear2 
perioarear2 multiarear2 spreadarear raoprogarear] 


Te Graphical Comparison of Model-Scale PSDs 


Plots of Model-Scale PSDs from all five methods for both the Wave Ht and 
TACS Roll are generated for visual comparison. From these comparison plots, it was 
determined that previously discussed “ptavg” selection of 50 was a very good 
compromise between display resolution and accuracy when compared to the other four 


methods of generating PSDs. (see figs. 7-10) 


oy 


%$ Graphical Model-Scale Wave Ht PSDs Comparison----- 
All 5 Methods 


figure (1) 

plot (wwelch, hwelchscaled) @ Welch Method 
hola son 

plot (wperio, hperioscaled) % Perio Method 
plot (wmulti, hmultiscaled) $ Multi Method 
DLOE( BURG Le) a runo 7 (223) $ Spreadsheet 
DICE tRUOwesGl > 47) 2RU67psa(s7 4). ® Rao Program 


title(’RO67CHO1 Model-Scale Wave Ht PSDs Comparison’) 
xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (dimesionless) ’) 

axis([0 6 0 .015]) 

Gira 


$ Graphical Model-Scale Wave Ht PSDs Comparison------- 
Spreadsheet vs. RAO Program only 


figure (2) 

Ot On) (eg ja tino 7 % Spreadsheet 
holdson 

DICE(ROG/pSa(27 47), R00 7psa ls; 179 % Rao Program 


title(’RO67CHO1 Modei-Scale Wave Ht PSDs Comparison’) 
xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (dimesionless) ’) 

axis({C 6 00065 ]) 

Gnid 


% Graphical Model-Scale TACS Roll PSDs Comparison----- 
All 5 Methods 

figure (3) 

pilot (wwelch, rwelchscaled) 

hOrGe on 

plot (wperio, rperioscaled) 

plot (wmulti,rmultiscaled) 

Dicwirumo7 (1,2), 7une7 (2,4) ) 

DlOECROG/DSda(>,47)7 R067 psa (+,.3,.); 

title(’RO67CHO3 Model-Scale TACS Roll PSDs 
Comparison’) 

xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (deg) ’) 

axcisi{[0-6 0.02) ) 

grid 


Welch Method 


de 


Perio Method 
Multi Method 
Spreadsheet 
Rao Program 


dP dO dO oO 
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& Graphical TACS Roll PSDs Comparison----------------- 
Spreadsheet vs. RAO Program only 


figure (4) 

joe cun6/ (:,2), Gpune7(:)4)) $ Spreadsheet 
Heweayon 

jumps OO 7/DSa(*,47), R06 /psdat{:; 3)) 6 Rao Program 


title(’RO67CH03 Model-Scale TACS Roll PSDs 
Comparison’) 

xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (deg) ’) 

eases i0 6 © .009]) 

fopanke! 


ABOARD 2 aia il YLT ENO Oe TENE LEI OLENA IEG Nt. saree rprnesmens remereear STAI oreronas gra macs ranges reenter et peer er tees ert SSE SI STR eae CSC aes oata ach tan ae aes 
M ; ‘ om "e j Fae ganel : 5 
: > o 


ROE7 CHO Model-Scais Wave Ht PSOs Comparison aa es ee ms 
_ Welch vs. Periodogram vs. Multi-taper vs. Spreadsheet vs. RAO Scan ae ie yf ; 





i' ois 


- | Note that the Spreadsheet and RAO Prograni 
-- .. 1 PSDs are much smoather while the others ate 
in their onginal, unsméothed state ' 


_ . PSD (dimesionless). 


: 





3 
. Omega ra 


vie : F eek, : ’ s xt. ae -_ St ES pts ‘ 
LEAS! MA tin Sat ER BIE sk oll L BRO Tee eh SteSetone 1: SIMRO Tree PR a Sat htc eth nc ahiace Merten ts all ARTE She Rae Gack eas Saige ak Giannis age Siena shia 


Figure 7. 


39 


+ 


bo 

ae 

ws 

4. 1 sielsislwiewieiwacesces sacs eocccccee 

ane : 

} 

s Jone we 

eo « ‘ 
“4. , 
"a 

Cxe ~ 


Fi 





= t+ is ar = » a Be ree a as wi eynte te qe Ba Sate es nee ut tee et awh 
r SY, a -_ - r oe 7 Mae, 1+ —  pEETCHOL ModetS Wo Ht PSDs Companson ©. “ yey ‘f ares seta i baer are age Mong ee 
' rs ee oe ee. el cale ove i bay ee TS. 
ay iol a OS ye cm > Spreadsheet Results vs. RAO Program Resutts "212" ne ine ee 
‘ és - I = 19 ud a se, o a & Libis exc! : “ae! 


erececeeeeeecaectoaaes 


eeereneeeeeeseousecne 


DR oes eemcnenenecseraasesetetoccseveceusaces 


) -- rk 


PSD (dimestonless 
i ; 


a“. 


ee Oe re ee ee 


ee ee er 


etait eT Oe ie tee el ete ete De el od ee ete eed Oe ee ee ee oe 


SSSeeeeeeesepesseeseeseeaeeepeseeson 





igure 8. 


Se twee Gere ee e2 eS Fer SSeS OF ereteGececerseseeraes sereeerereaes 


‘RAO Program PSDs 


‘Spreadsheet PSDs 


Se eee Pe eese See ees SOF eHFeeS SESS Sreseserseeecue 


eececeteoee 


eeeee ee eee seuseeesecen 








SOP OS SOSSHSHHSESSSSES SESH HOEEH SEE SEBS OPE eeZeeseene 


| 


rerTrer er eer? yt ee 
eer -eeee-eree 


SOSSHS SESS SSHSEHSSSHSHSSH SHS HOHE SeSeBESS 
. 


Seer erase eeeweeeSeeeeeeseseeseseeeeese 


ee ee ee oe ee 2 ee es 2 


e®eeeeeseeeeeee races eeeneececoesesceseseeeeeseeeeeesectese 


ee ee ee ee ee eee eee eee ee ee eee 


~ 


errr ereree-ehee tere “FS eeetehoewr--e 


steerer ee-ees 


i tak 
Bia FiO 
Ali 2 ad 





” 


. (Wer . . . ° ** 
, A J i ‘ me rie Nae eae re ‘ Oe ae ad * 
Tee cacgrens 3 Tha oat canara Teton wae! Pome KI re UN (a TPR alia ic chest TEESR: Seis. tS V Cr cctsecBalsh tn Mes iad cd Cenait bet Recc be sei BR PRE iy RITES IR Ra ae 


is 


« 
- 


Tar EST, | rem enrnee ~~ ROB7CHO3 Model-Scaié TACS Rol PSDs Comparison “= eae ee Sn ame ae Se an eae 








F ) e 
2 ° ' 
d ° ‘ 
z s s 
‘ ‘ 
s s 
a7 s 
° . 
ob ce ccecccccccccoccccsocced 
s 
e s 
4 es s 
e ° 
1 e ° 
° ® 
e s 
s s 
s ' 
s ‘ 
OSES OOOO OT OSOs Oe BESS ee aneereeeeeeee 
° ° 
e ' 
es s 
e e 
vs ° ° 
e . 
e ' 
s s 
e s 
& ones ecccccccca-secseeced 
s . 
° . 
& ° e 
s ' 
rs s s 
s . 
° e 
° ° 
° 
cad s 
e s 
A s ' 
* s ° 
° ° 
s ‘ 
e s 
e e 
s es 
s s 
Saeed OOOO oe Ce eee eee ee Oot cobs ceeeceseees eeee cesses eseseeeeeeees, 
e ° 
an, ; : 
. ° 
s ' 
e ' 
s e 
° ° 
. ° 
a. : : 
Oe ee ea 
e s 
s s 
e e 
° ° 
e e 
e s 
e. e 
e s 
s e 
a Ooo Oe 08 OOS HOH OO 86S OOM SSO SOOO OSSOSSESS OS Sere MESES OE 
+. . ° 
7 . e 
. ry 
< ry ° 
3 e ) 
e s 
e o7 
. ° 
° ° 
OO Fe oe Or Oe Seer sees ee coteesr cess Benen er Se eescadagacse 
« 4 
e ' 
e ' 
e e 
. ° 
e e 
s ' 
° ’ 
. e 
‘ e e 
‘ ° 
, 4 
4 
‘ 
e 
4 
° 
4 
‘ 4 
a) 
Biv eos nyse yam de ele Soy sti wn ween aaBisano 2 cE Enisteranbenaetin © anak Low Kine d te Petadmeatndiere tian tll hin ry drent eh 


Freure' 9? 


ee. apa es ” Welchvs, Perlodogram vs: Mult-taper vs: ‘Spreadsheet ys. RAO Program ~ 4 /Beg, ae = a 











rs oe 


hh ne HE Ree 


SOP eeccesocecteccoeeesoccosetlece sees cece e@esecees ease Oheteoetoeceacesets see 88en 
. 


peccccccegtgecee=- 


OOO O88 Se OOO SOOO ESSE OS OS SHOE OES SHOE OS SEES EH HESS SS SS HS SESH Sle SESS SS SeseLeeeeeus 


Pet ecccegtetocsecse-Seee-seee 
tem geeceggert eereee 





serer-ccoecpe 


Dens ccasassiccsteesaeccocs{ccccdacesccsesscescsccsch coeeoeeescuseeaeoseres 


Note that the Spreadsheet and RAO Program 
PSDs are. much smoother while the others are 
in their original, unsmoothed state 


Oddcb ues concsceceescasceS cee scoecceeee ce eeceeccechkiccécleceecceeceeee esa es 


eevee cetecscceg-gtgeg-e-secebessceecasy 


Or ee ee 


eere-esege 


e 
e 
. 

° 


© 00 we ee eee eee eee ce £80008 ose se oses esses esse ecess Sh eee Sesser eser eectea=eee 


SOSH ooetetSSSSFOSZe*ESHFSSESES MeSegeeege 


epeeeertere 


~ Omega (rad/s) . ae fio Oixs ~ ae 3 - we , ' 


Ayr \ « 


Po po oe A RN I, De = Pe UR Gn OO WATE, OY Sup Pee ne TRV Wee UCU ov «NY 


40 


ree coremee 27 


| “posroHs ModelScele TACS Roll PSDs ‘Comparison 
Spreadsheet ve” RAG Program Only ee Oe 


neers Spreadsheet Rell. RAQs.......... 
RAO Program Rall RAGs 


eeweeseeescesesensensheccceseeeepesetecanamatans Pannecsveseenceescbtanecseecenseeesseese en bone aese ne eeeseeneeneae 


eee ee ee 2 


‘eee om ad ae 
—— BN e. 


Ss => 432 


BS 
« 


CeO OR eee OS OO SEF Sees sesenFeseseseceseeeeseevease SPP sess sees Hfoce deca Fe SSF SPSS SF SSHSSF HHS SHSSSHTESSHSSSETOSHSESOHH BES EBeeseseesea 


7 
Or 


Ore 





psp taal 


Cee mere ee eee PHBH os See EFF OSSeeesesF®eseseeese ease esesenennasee2eeen 


rere ee err ee ee eee ee eee ee eT CTT TET CTT Tee Tee 


coow Vi gw gape ete e se eee sees Fe esses eeseseeseeeGeneeseeeese Seeteocee 


wmoeceeracgrar 





a 


ey ona ’ 74 “| ae “ tutte’ mel A . . , ¢ ‘ Olas. a t Ld af wf . ° oa 4 oe 
eae a. ee oe Oe. ee v - i ae” a Se Le eg 1 a | ee Mee 
Dhara! Decinast Sree int i Lcreechalies: Zee ine Re Dro SoA EA a Ee I OED PIO Dy Pad PEP IOP: PN ae a ame. sate Pi Ob ts wt Nee ban AP teen aan Kath, on Mle cht SG eaten A ate AB, de eS Sek, Gece salah rei mele alts cee eat ate: Util tle Pag isdn % 


Figure 10. 


8. Compute Model-Scale RAOs 
Just as it is performed in the RAO Program, Model-Scale RAOs are computed for 
the Welch, Periodogram and Multi-taper methods by taking the square of the scaled 


TACS Roll PSDs divided by the scaled Wave Ht PSDs. 


$ Compute Model-Scale TACS Roll RAOs------ eo---------- 
Ri = rwelchscaled./hwelchscaled; 

welchmrao =) SGuE(R )- 

R2 = rperioscaled. /hperioscaled; 

periomrao = SGEE(RZ); 

R3 = rmultiscaled./hmultiscaled; 

multimrao = SOrE(RS) 


4] 


Graphical Comparison of Model-Scale RAOs 


2. 


Plots of Model-Scale TACS Roll RAOs from all five methods are generated for 


d they provide very satisfactory results. (see figs. 11-12) 


| comparison an 


visua 





AEE AMEN TS CSREES Of CUARY CREOLE RR oN CRED SP COC heat Sag MANOR On met AARDT, BR ES tO ERMAN 
Rie SOR a es SRS RAE AR Bas BAERS Se 

hb os te x : Ey! 

JX) z ” gis 









Sow OSS ESOS EGE SE SEEEES EMSS SESESEER SETH HOSE EESESSSSSEE BP eeeseesestesesessesgewesessses 


¥ 
“29 gee AWE "Caen Ey ‘cats Mae, Wr RR ie SI RS 
vy ee re -t = ar p> a e . at a x eit a wy gh x 5 at? Sve 4 if a robs % 
a EN eee Be 5 Ug U0 BORN BRR OE Se ERR RE RUE Bie BEY a, ry Bhante 
Pa @ on oft AGT Ehas Rc ee AN pas PE ne Be eS ee ie, soli > aS eS x} 
Se eS ta a acetate ag SOR eS TRE SEE oh eae ERS 
aa * 
e ¢ ry 4 ~ es 
‘ e Py ——h et 
r ‘ ‘ ‘ < 
: : ae 
U ‘ ° ‘ Sha 
: : : : —{ ie 
: : 4 
: e s 1 aie, 
: ee 
: ‘ e ‘ ‘ a3 
: : : : : P|: 
5 a ® e ‘ mo 5 ie 
Red : : : : = ; 
ec Og ‘ e ® ‘ 4 be 
San eee Pees Sent en erent ore r eres fn newelsice deeccccccerPeccccnnes dececccecce becceee A 
Fe ° ® ‘ 573 
ae ‘ ° ‘ e a ¥ 
NS ® ° ‘ e eI 
; ets : : : 5 ; rs 
peg 7 e ‘ ° , .. 
— ‘ ‘ ‘ J 
eo : : ° 3 
rie> o . r) r . e Ra 
s ’ 
A = Wl e * aoe 
aN 8 cy ® a 
, 
hei : : - 
E e e 
* . s s 
20: : : 
Rae ° : 
io : : 
“a4. 


ae 
‘3 
83 
& 


te 


anne 


‘~~, 


POPP eo ea ee ewes eee FSSSSSSSES ASSESSES SESE HP SESE SSESS|ESSOS St ewerwisseseeseegesesesene 


a ad 
ss 


b-baper vs 


OOO Oe OO Oe 08 OOOO FOO OOS OO OS OOS HOO ESET OOS ES O85S OSES FOS O42 55 28 288s 228888 S88 


2808 08 O48 O06 6088 O88 298 888 S58 COK 08 28S £8S8 88284 Oe eee oS Beebe ee es CE SeeSeseeseesess 


eeee cosdooqesezcensccscscedesssctgecsese 


OO | 


\ 
~ Omega (radi) sor 


ut ae Sy 
, “an eho Eee (AES TRY Wn gaw re 


Sf 
oy) 
” 


Mu 





“« 
» A 


c, 4 
ve Arig Lota: take 
> 2 ie 





ram vs; 


s he 
Ey Ay 
sé 


ne nccecenencncerenc ence ndencececcncencnccneescenes 


gram 


Sa 





“ROB7CHO3 ModE Stalé TACS ROM 
9 











s 
. 
8 
s 
s 
e 
s 
g @ : é 
CT LAY) ® 8 
ewecrece= aya 
Oo ooee- Bild eee eer ect e tec c eg cee age Bore Seem cence cc cc ccc cn ges coscccoeqeccccccccs NYY 
"tc s e i] ry ? 
a iS o as 
‘a. £ ; : S “ye, 
Qa ) ra) ° ® e ae ci 
2p ‘Oo M : : Si of 
2 i 2 : : ° Soa 
it~ £& ® H : : Na 
P= 1 @ @ ‘ e ‘ ; tet 
2 ‘Os % ' : : eony 
Cc : . ® ses 4 
£ 1BEo : : ; : wets: 
: ac : : ; : ae ere 
r nn Nee +e 
ri Ne ® & ' PY ® ® ° 2 Oa 4 
. = 8 8 ; : ; : ; 7 
3 wn ‘ ® e . ——— es 
pe dk Preece eee IE) -O f..1.-..----- Joe ccccccccbecesccees Beceeee eens Mee eeccceee Pecececeee a _ ms 
= ¢ VEY ‘ ° ® ‘ a 4 t. 
Le sO we c : e 4 Ae 5% 
“ = s e 8 eo. Pac 
: Pa ; ‘Oc = ‘ e e ° e 25 & 
i Hon : : H H ® a 
Ege i . =~ AY . ® e ‘ Py 4463 
 - Marsh E Cc . ——— eee ~% ’ 
wt ‘GS On ry e ry) r ry = sy 7 k 
-. ie, DS : : : = SS = > ' 
; ‘ i= * ‘ ry === + ae 
Fi Sg. aS) \y) Qo 6 : . 5 ae 
‘eo ‘Sw = : ers 
ee ae) ® a 
i 3 Ke) < = Py coi 
et 2 fi 
Baad, to ; ns 
¥ =e 7 ef 3 
A om ‘ . go 3 . 
fe 4 . a ; y ae ae 
a - ee Se «84 
2 pet witb Ne as “ 
ee BS Rio f f | 
5 <S 2 5 3 = | “ ag aly - Sabo ae be Rane $ ef 
SRE I BA RERES F BADY RE NOOO ee eg. SECT EWAN, Late a Sad ets med 


Figure 11. 


42 





* “ a 


© roszcte9 ModelScale TACS Roll RAOs Comparison’ ~ é 
Spreadsheet vs. ‘RAO Program Only Ag 





wees ase seeeeeesceeebew es eee eee sense eteeesetbaneececccesceceeeeenececeebiaweceecsc ees eceeeuseanmmtboenseneeeseeeseeeneae 


eee er eeeeee 


Spreadsheet Roll RAOs 


RAO Program Roll RAOs 


eeer eset ote eter eesHeeeeesea 








Sea detiate bleed AD Ct Tey See eT) 
Se OS OSEMBOPP MESSE SEEHEEBEEES ES 


Sete eee oF SSeF OSE He gheroeee SS Ser SeSSSeeseerebeeaeeW=FF OFF S888 SOB eheaBeeesSSSEGeneeseese 


 ? fe eee eee eee eee eee eee Teer rT ee Te TT TTT Te 
Pe ee ee ee eee rer Pierre TTT TTT TTT TTT Tee eee 


Ce ee ee ee 


SBoocnccoogscocoos ese, 


as, Ke Sie, 
oe Eas 

oe a > ae, 
fey: oe “ss, % 5 wh es 
‘<3 He 2, 





‘ <a Fe 1 £& 4 Sea P wie “yw. - % eS 26 awe. ae hs 
55 BitS aFetA osan SaccteAnina obec Sass Namen eke fe ASS er RAS <2 Saeselstenn ee PAE AAC OD a A MET SEM OR OU Ae ST TT RR 


Fis gure 12. 


10. Scale to Full-Scale RAOs and Omegas 
In order to compare Full-Scale RAOS and Omegas, the Model-Scale RAOs and 


Omegas are scaled accordingly just as they are done in the RAO Program. 


@ Scale to Full-Scale RAOs and Omegas----------- ------ 
welchfrao = Wwelehmnag-w!2 724.175. 
periofrao = penlonmnac- 12772454 75° 
multifrao = Mule Mmao 7/7 24 no. 


wwelchfull = wwelch/saqrct(724.175) ; 


woperiofull = wperio/sqgrt (24.175), 
wmultifull = WMUlGl, SQ (74.475) ¢ 
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Jt Graphical Comparison of Full-Scale RAOs 
Finally, the Full-Scale RAOs and Omegas for the Welch, Periodogram, Multi- 


taper and RAO Program methods are plotted for visual comparison. (see figs. 13-14) 
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Iv. CONCLUSIONS & RECOMMENDATIONS 


A. CONCLUSIONS 


Due to the validation comparisons performed on run 67 above, the RAO Program 
performed quite admirably and accurately. The RAO Programs ability to accurately 
produce Model-Scale PSDs, Model-Scale RAOs and ultimately Full-Scale RAOs for 
each channel, along with their associated Omegas, attests to the programs performance. 
In addition, a short looping program, “raoloop.m’, enabled up to seven data runs to be 


evaluated sequentially. 


$ Short Program Title: "raoloop.m"-------------------- 
eum [1921 23 ee oe 1 ee $ DTMB Run Numbers 
LOR sZea= 1 > / 

rao(a(1,zZz)); 
end 


The seven run looping limit was imposed due to local hardware configuration 
limitations. With each data run producing approximately 12.5 megabytes of stored 
figures and matrices, a 100 megabyte zip disc could safely store up to seven data runs. 


With a read/write CD drive, the looping capability would have been much greater. 


Upon completion of processing the Configuration I runs, the composite matrices 


(Model-Scale PSDs, Model-Scale RAOs and Full-Scale RAOs) for each run were 
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assembled all on one compact disc as a convenient reference from which any two 


columns of results may be plotted against each other. 


B. RECOMMENDATIONS 


An opportunity for a logical continuation of this thesis will soon be available. 
Trial runs using the actual, Full-Scale ship, the USNS Grand Canyon State (ex-SS 
President Polk) (T-ACS 3), were conducted September 9th-16th off the coast of Southern 
California near Camp Pendleton. Data recorded from the Full-Scale trial runs could be 
processed in a manner similar to the RAO Program and used as a comparison tool to 
either disprove or further prove the validity of results from the RAO Program. If RAO 
results from both the at-sea Full-Scale trial runs and the RAO Program model-basin 
evaluations are in agreement, a multitude of combinations can then be evaluated to 


predict ship motions in an irregular sea by using the principle of linear superposition. 
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ua 


G25 a2 54.5 





paar 
w> 
I 


b ie) 
wy 
| 

= 

oe 

n 

Q 

Ss 

~) 


me ae) 
wf ~) 
I 








444 AUG = 


14- 
a] 


14- 
448 Aug- fTASCA7 447,449 190 iI 180 
97 


270 40 0 54.5 


14- 


TASCA7 449,451 197 





453,455 JL Rew [16 180 


180 


=f ffs 

=F FET 
315 0 b fs 

a 
=f Fe 


o5 20 


Aug- ITASCA7 455);457 194 
97 

14- 

Aug- JTASCA7 457 9S 
97 

15- 

Aug- gTASCA7 459,461 19641 180 
97 


456 
458 
462 1971 270 
15- 
aug- |TASCA7 463,465 4198]T 180 
97 ) 
15- 
aug- |TASCA7 =«-*[465,467 = 199 fT 270 
97 
467,469  |200|r 180 
469,471 {201 Ef 


180 





S) NS) rss rs ww wr ) 
+ w 





Aug- gTASCA7 F 


vat 
es 
6 9) 
0 DP 
ii 
QQ 1 
! 
ry 
> 
p 
~) 






470 Ss 
474 TASCAT 473,475 203 iI 200 5 0 0 25 






NREL Data Lost. 


MEO Dh Pp 

in Mae gett 

i iQ | 
\ 


a5 


+ 


4) 


176 so-foascar fons. oaks fe 
13- 

478 Aug- JTASCA7 477,479 205 {I 180 Eyes) 
97 


TASCA7 479,481 as oi5 


‘Oo 
4 


INRL Data Lost. 





bd 
oO 
re, 
oO Db PR 
Ic ul 
QQ | 
\ 


15- | 
482 ug- [TASCA7 481,483 207] 270 14+ 4315 25 
97 | 
TS 
484 Aug- |TASCA7 483, 
97 






a 
\O id 
~~] ul 

1 


Data Lost. 










ITASCAT 489,491 placa Wk 
. INRL Format 2. 
TASCAT 495,497 21451 4180 815 29.25 ly 25 
, INRL Format 2. 
TASCA7 499,501 216 {I 180 els: 29329 12 25 
: NRL Format 2. 







22541 






5g 


THIS PAGE INTENTIONALLY LEFT BLANK 


60 


APPENDIX B. TEST CONFIGURATION II 


INRL [Test ave ‘Séa Boom 
onfig Head State Slew 
# ing Angle 
15- | 
hug- ITACSA7 218 4Ii 
a7 al 
15- 
506 Aug- IJTACSAT 50D, 207 219 4Ii Z/0 ae B15 5 iz 25 
“i | 
S709 220 }Ii 


DNSWC 
un # 






\s) 
i) 
ct 
0 





wm 
O 
rey 


fo 

OO 

Oo 

vam, 

tJ 

pay 

ut 

WM 

fe * 
NO 

rw) m 
ui 





508 





a 


fe 
ut 
t 


Bi 2 


‘oO b 
Ie 
(Q 

{ 


Ut 
7) 
oO P 
Je wm 
Qa 1 
{ 


223 


515 $14,516 {<°- |r 90 2B 270 NRL Load Data 


Not Locked 


Paite= 
S17 Aug - ITACSA7 56,5 L8 224 ]4Ti 3s 270 5 a2 25 
97 | NRL Format 2. 





OP 
IE Oo 
Q 1 

{ 

fe. 

oe 

Q 

Be 

i) 

a 


521 Nug- ITACSA7 


W 
=) 
iO 

oO PF 

JIE a 
Q | 
\ 











































523 $522,524 erun to get 
good load data. 
526 TACSA7 529 INCD jIi Load Control 
| Configuration 
S27 ASCA8 525 NCD ITi Load Control 
Configuration 
528 ASCA8 S25 NCD {Ii 90 Road Control 
Configuration 
ransition 
529 TASCA8 525 NCD|Ii Gs aaa 
Configuration 2 
| EOS 
530 ASCA8 525 INCD fIi Load Control 
Configuration 
Dynamic Load 
a TASCA8 525 INCD ITi Control 
| Demonstration 
32 TASCA8 $25 NCD {Ii Load Control 
Configuration 4 
Dynamic Load 
555 \ Demo No Rider 


Block 
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APPENDIX C: TEST CONFIGURATION III 


R Test ave 
ot poof psa 
F fh State 
Directory # ing 
re | 
Soil Aug- JTASCAY 536,550 bia 
| 97 
17- 
539 Aug- ITASCAY 538,540 22 7aiai 
97 | 
540,542 2284111 
543 Aug- ITASCAY 542,544 229 4iii c’ 
97 
17- | | 
545 Aug- {ITASCA9 544,546 230 tii 3s 
547 
549 Aug- ITASCA9 548,590 1S 271 1 1 p 
TASCA9 550,292 2350111 bo 








w 
i 







i 
‘s E 


WwW 


541 

















Sisal 
= pascas fsa. ss 4 Ae pf 
Bes TASCAY 235 43:3 
iv pascas_fise.se peop i 
e. pascas fis, se 
Sa TASCAY 560,562 
1563 paseas fea. se 
Bes TASCAS 564-566 

97 

17- 
568 Aug - 

97 | 

en. 

Aug- JTASCAY S71, 573 

97 


wm 
~) 
thd 


74 fiv- frascas_s73,575 pasfiar__[s0_p [0 [s3 fa feo] 
oye: 





ug- 
97 
17- 
576 hug- {TASCA BI 52577 2444111 3s 2 Sas 13 
97 
17- 
578 AUg-ITASCAY 5777579 245 fiii 3s a5. 3 13 
97 
17- 
580 hug- JTASCAY 
97 
Sea 
Jo3 hug- ITASCA 582 
97 
17- 
584 Aug~- ITASCA9 582 
o7 | 
17- 
565 hug- {TASCA SoZ 25040 1 2 
97 
17- 
586 hug- ITASCA9 Buc ii] 
97 














240secs Static 
Load Control 
Demo 


240secs Static 
Loaa@ Control 
Demo 


240secs Static 
Loaa@ Control 
Demo 


200secs Static 
Load Control 
Demo 
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APPENDIX D: 


Channel! # 


MASTER CHANNEL LIST: CONFIGURATION I 


DTMB Runs 1-503 and NIL Runs 1-217 


Channel Title 
Wave Ht Bow 
Ssonix Sonic 
TACS Roll 
TACS Pitch 
TACS Roll Rt 
TACS Pitchkrt 
TACS Yaw Rt 
TACS CG ZAcc 
TACS CG YAcc 
TACS CG XAcc 
TACS CT ZAcc 
TACS CT YAcc 
Cnir Roll 
Cntr Pitch 
Cntr CG ZAcc 
Cntr CG YAcc 
Cntr CG XAcc 
TACS Cnir Rel X 
TACS Cntr Rel Y 
TACS Cntr Rel Z 
Lhtr Roll 
Lhtr Pitch 
LAtr Roll Rt 
Lhtr Pitchkt 
Lhtr CG ZAcc 
Lhtr CG YAcc 
Lhtr CG XAcc 
TACS Lhtr Rel X 


TACS Lhtr Rel Y 
TACS Lhitr Rel Z 
Load Mo Orthognl 
Load Mo Parallel 
Boomflip-HorizAcc 
Boomlip-Vert Acc 
Lghtr4-Strn VAcc 
Lghtr4-Bow VAcc 
Lghtr3-Bow VAcc 
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Units 
inch 
inch 
degree 
degree 
deg/sec 
deg/sec 
deg/sec 
GS 
GS 
GS 
CS 
SC 
degree 
degree 
CS 
CS 
CS 
inch 
inch 
inch 
degree 
degree 
deg/sec 
deg/sec 


Final File Title 


r###CHO!] 
r###CHO2 
r##ECHO3 
r#E#CHOA 
r###CNOS 
rHEECHO6 
r#H#CNO7 
r###CHO8 
rHHECHOS 
r###CH 10 
rHEECH 11 
rHE#ECH12 
r###CH 13 
r#H#H#CH 14 
r#E#CHI5 
r#HECNI1O 
rHEHCH | 7 
r###CH18 
r##ECHI19 
r##ECN20 
rHEECH2 | 
r##H#CH22 
r#HE#ECN23 
rHHECH24 
r#E#CN25 
rHF#CN26 
r##ECH27 
r#HE#CH28 


r#EECN29 
rHEHCHSO 
rHHEChHS | 
rHEECH32 
rHEECH33 
r###CH34 
rHHHCH3S 
rHHECN3O 
rHHECHS/ 


Lghtr2-Bow VAcc 
Lghtrl-Bow VAcc 
Lghtr2-SMdshpVAc 
Lghtr2-PMdShpVAc 
Lghtr2-CLMdSpVAc 
Lghtr1-CLMdSpVAc 
Lghtr4-StrnLngAc 
Lghtr4-PStrnivAc 
Lghirl-PBowlvAcc 


66 


YOODODODODOOO® 


r###CH38 
rEEECH39 
r###CNH4O 
r###CHA) 
r#HECHA2 
r###CHAS 
r###CHAA 
r###CH4AD 
r###CH46 


APPENDIX E. MASTER CHANNEL LIST: CONFIGURATION II 


DIMB Runs 504-533 and NRL Runs 218-226 


Channel # Channel Title Units Final File Title 
0] Wave Ht Bow inch r#HECHO! 
Q2 sonix Sonic inch rHEECHO2 
03 TACS Roll degree rHEECHOS 
O04 TACS Pitch degree rHEECNOA 
O5 TACS Roll Rt deg/sec rHEECNOS 
06 TACS Pitchht deg/sec rHEECHOS6 
07 TACS Yaw Rt deg/sec (HEECNO/ 
08 TACS CG ZAcc S r###CNOS8 
O09 TACS CG YAcc S r#H#ECNOO 
10 TACS CG XAcc GS r#EECNIO 
1] TACS CT ZAcc GS (HHECHN 1 1 
12 TACS CT YAcc S rHHECHI2 
le Ltr Roll degree rHHECHI3 
14 Lhtr Pitch degree rHHECNIA 
is Lntr Roll et deg/sec rHHECHIS 
16 Latr Pitchet deg/sec rHHECNIO 
17 Lhtr CG ZAcc GS rHHECHNI 7 
18 Lhtr CG YAcc GS rHHECHI18 
19 Lhir CG XAcc S rHEECHI9O 
20 TACS LhtRelX inch r###CN20 
Z| TACS LhtRelY inch r#HECN2 | 
22 TACS LhiRelZ inch (HEECN22 
23 Load Mo Orthogni inch rHAECN23 
24 Load Mo Parallel inch rHEECH24 
25 Boomflip-HorizAcc S rHEECN25 
26 Boomflip-Ver Acc S rHEECN26 
7 Lghtr4-Strn VAcc GS (HEECH2/ 
28 Lghitr4-Bow VAcc S r#HECN28 
29 Lghir3-Bow VAcc S rHEECN29 
30 Lghtr2-Bow VAcc GS r#EFCHSO 
$1 Lghirl-Bow VAcc GC rHHECNS | 
32 Lghtr2-SMdShpVAc S rHHECH32 
33 Lghtr2-PMdShpVAc GS (HEE CH33 
34 Lghtr2-CLMdSpVAc S rHHECNS4 
aS) Lghtrl-CLMdSpVAc S rHHECNS3S 
36 Lghtr4-StrnLngAc GS rHEECNH36 
2S), Lghtr4-PStrnTAcc GS rHEECH3S/ 
38 Lghtr1-PBowTvAcc S rHHECH38 
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APPENDIX F. MASTER CHANNEL LIST: CONFIGURATION III 


DTMB Runs 534-587 and NRL Runs 227-25) 


Channel # Channel Title Units Final File Title 
0] Wave Ht Bow inch rH#HCHO! 
O02 SONIX SONIC Inch (HEHCNO2 
03 TACS Roll degree (HEECNOS 
O04 TACS Pitch degree (HEECHOA 
05 TACS Roll Rt Gea sce r###HCHOS 
06 TACS Pitchkt deg/sec r#HECNOO 
Q7 TACS Yaw Rt deg/sec (H#ECNO7 
08 TACS CG ZAcc SG r###CHO8 
09 TACS CG YAcc GS r###CNOS 
10 TACS CG XAcc GS r#EHCHIO 
i TACS CT ZAcc SG rHEECN |] 
12 TACS CT YAcc SG rHEECH 12 
le DD 963 Roll degree rHHECNI3 
14 DD 963 Pitch degree r##HCHI4A 
is DD 963 CG ZAcc S rHEHCHIS 
16 DD 963 CG YAcc S r#HHCH IO 
17 DD 963 CG XAcc S r##HCH 17 
18 TACS1 DD963 RelxX inch rHEECHI8 
19 TACS1] DD963 RelY Inch r##ECNIO 
20 TACS1 DD963 RelZ Inch r###CN20 
Z| Load Mo Orthognl Inch rH#ECN21 
22 Load Mo Parallel Inch rHEHCN22 
23 Boomlip-HorizAcc S (HEECH23 
24 Boomlip-Verf Acc S rHEECN24 
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APPENDIX G. RAO PROGRAM 


PMC tLOMm raotrun ) 


fold 


Program Title: rao.m------------------------------------- 
Referred to as RAO Program 


feed 


oe 


Must enter one of the following to initiate this program: 


feed 


fo) 


rao (##) for DTMB run numbers <100 
rao (###) for DTMB run numbers >100 


oe 


% Initial Set-up------------------------------------------- 
sread = char@s 
swrite = (ena Gabe 


format short e 


ptavg = 50; 
window = S008 
lambda =e ee Les 
freg = 2c 


% Establish Proper String Configuration for Data Run------- 


(o 


Vie saibigmeec— S/S, 


Sk =) Chaat” RO). 
$2 —Beelsisue (4) 
else 
Sul = IGChart rR’); 
end 
srun = numMzses (rium, 
filename =" streatlel,srum) 


Perri -<=eoo 


folder = “Streaie(sZ,srun)- 
else 
folder = Sri 


end 


aA 


% Determine Number of Data Channels to be Processed-------- 
ge G Urals — as 0) 


channels = 46; 
elseif run <= 533 

channels = 38; 
else 

channels = 24; 


end 


% Establish Filepath and Load the Raw, Time-based Data----- 


filepath = {[sread, filename] ; 
load(filepath) 

szdata = s1ize(data) ; 
channels =. 'SZ¢data (Zz); 

n = szdata(l); 


%$ Compute both Model-Scale and Full-Scale Omega Ranges----- 
momega = 2*pi*freq/n* (0:window-1) ; 
fomega = momega/sqrt (lambda) ; 


6 Full-Scale Scalers, PSD Units and RAO Units Library------ 
fsl = 


pul =f (aan 

rae oil = ’(dimensionless) ’; 
fs2 = 12/lanbda; 

pu2 = ‘(deg)’; 

ru2 =~ * (deg/in) ; 

fs3 = j12/lambda~ (3772): 
pu3 = '’(deg/s)'; 

rus = '({(dég/in s)’; 

fs4 = 12/damoda 2 ; 

pu4 = ‘'(deg/s%*2)’; 

ru4 =" (deg/ian ss" 2) 

% Looping/Processing of Data Channels---------------------- 
for ch = 1:channels 


if channels==46 


ame o—— 
fscaler = 
psdunits = pul; 
racunits -= rul; 


F2 


chtitle 


elseif ch== 


fscaler 


psdunits 
raounits 
chtitle 
elseif ch== 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==4 
fscaler 
psdunits 
raounits 
Cherelle 
elseif ch==5 
fscaler 
psdunits 
raounits 
chtitle 
elserr ca—— 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==7 
fscaler 
psdunits 
raounits 
Gheaede 
elseif ch== 
fscaler 
psdunits 
raounits 
chtitle 
Srseickca=— 
fscaler 
psdunits 
raounits 
Ghtitie 
elseif ch==10 
fscaler 
psdunits 
raounits 
chtitle 


‘Wave Ht Bow’ ; 


=r Si 
Pu, 
ru 
SSOnixesoomilc 


fs2; 
pu2; 
cue: 
SACS RO li 
iso2- 
pu2; 
1g, 
TRACS bite. = 
Ess 
pus: 
IGACLS) 
ABE CIS, aieeyl Gl Sie, 
Eos 
puss 
Euler 
"TACS 


fs3; 
Pus: 
ruse 
"TACS Yaw Rt’; 
£s4; 
pu4; 
ray 
TRESS CG LACS; 
fs4; 
pu4; 
rudd: 
“TAGS CG: Yee; 
fs4; 
pu4; 
ru4- 


“ERCS CG xAec’- 


i. 


Plten Re = 


elseif ch==11 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==12 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==13 
fscaler 
psdunits 
raounits 
Gntertile 
elseif ch==14 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==15 
fscaler 
psdunits 
raounits 
chtitle 
elseif -ch==16 
fscaler 
psdunits 
raounits 
cntaele 
elseif ch==17 
fscaler 
psdunits 
raounits 
chtitle 
elseit “ch==18 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==19 
fscaler 
psdunits 
raounits 
chtitle 


Esa 
pu4 ; 
nag Os a 
‘TACS 


fs4; 
pu4 ; 
ba Dil. 
‘TACS 


fs2; 
pud; 
ie 
“Crler 


fs2; 
pu2 ; 
sal Baye 
eB: 911 Bn a 


Loa; 
pu4; 
mua: 
Cree 


fs4; 
pu4; 
nud: 
sig i 


fs4; 
pu4; 
mud: 
fe Bie Hwa 


fsl; 
ober 
nag We 
"TACS 


Esl; 
Duel; 
Tal; 
‘TACS 


CEIZACe > 


Ch tACe.s; 


Roll <- 


Pitch’ = 


CG ZAGcc > 


CG YAcc’; 


CG XAcc’; 


Cntr Rel X’:; 


Cher 


Rel Y’; 
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elseif ch==20 
fscaler 
psdunits 
raounits 
chtitle 
ellicier ba ehn=—2u 
fscaler 
psdunits 
raounits 
chtitle 
elserf  ch==22 
fscaler 
psdunits 
raounits 
chtitle 
erserh.ci=—2 3 
fscaler 
psdunits 
raounits 
chtitle 
elseith ch—-—24 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==25 
fscaler 
psdunits 
raounits 
chtitle 
Sls eat peh==z 6 
fscaler 
psdunits 
raounits 
chtzetle 
Slee schi——2/ 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==28 
EFscaler 
psdunits 
raounits 
chtitle 


fsi1; 
pee, 
ue 
TAGES 


ESZs 
IDOL, 
Guz; 
Sianess 


See 
[OO - 
ruZ 
IU IgiEse 


Ess; 
Dus 
Eur 
* eae 


Eso 
Dus. 
Tatlin: 
’ Totaiteas 


fsa; 
pu4 ; 
Tae 
’ Deis 


fs4; 
pu4; 
ru4 ; 
ge oleae 


fs4; 
pu4; 
rude 
eo lepe 


fsl1; 

pues, 
rul; 
‘TACS 


Ghee | yo: 


Reis 


Paeeia 


Piten 


CG YAcc’; 


CG XAcC’ ; 


iiiiete kel 


/ 


Roll Rt’; 


RCo 


CG ZAcc’; 
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elseif ch==29 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==30 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==31 
fscaler 
psdunits 
raounits 
Chettke 
elseif ch==32 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==33 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==34 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==35 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==36 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==37 
fscaler 
psdunits 
raounits 
chtitle 


fed 
pie, 
iol ibe 
TAGS aie Re | Yo 


ks; 
pu: 
Gu; 
TACs LHNtr Rel 27s 


feds 
Duds 
ieee 
‘Lead MoyGrhnognia- 


Esl: 
pul; 
aio (Ale. 
‘Load Mo Parallel’; 


fs oe, 
pu4; 
ru4; 
‘Boom: Tip-HorizAce 


boa: 
pu4; 
ru4; 
‘Boom Tip-VertAcc’ ; 


£s4; 
pu4; 
ru: 
‘Lohtr4-Strn VAcc’ ; 


fs4; 
pu4; 
ru4; 
‘Lohtr4-Bow VAcc’; 


fs4; 
pu4 ; 
rua: 
‘Lgohtr3-Bow VAcc’; 
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elseif ch==38 
fscaler 
psdunits 
raounits 
echtwt le 
elseif ch==39 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==40 
fscaler 
psdunits 
raounits 
ehrart le 
elseif ch==41 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==42 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==43 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==44 
fscaler 
psdunits 
raounits 
Cheat Je 
elseif ch==45 
fscaler 
psdunits 
raounits 
ehe wile 
else 
fscaler 
psdunits 
raounits 
chtitle 


f£s4; 
pu4; 
Weve 


‘Lgohtr2-Bow VAcc’ ; 


Esa 
pu4; 
TaN hs 


‘Lgohtr1l-Bow VAcc’; 


fs4; 
pu4; 
nu4- 


beGhtmaZ2—SMaAShp vee 


fs4; 
pu4; 
jeV YA, & 


‘Lgohtr2-PMdShp VAcc’ ; 


jel 
pu4; 
ia 


‘Lghtr2-CLMdSp VAcc’' ; 


fsa; 
pu4; 
ru4; 


Egan CiMasa Vece = 


fs4; 
pu4; 
ua. 


“Lentr4-Strnimg sace: 


fs4; 
pu4; 
rie: 


EGHNEra = PStrin EVACG: ; 


fs4: 
pu4; 
rua; 


‘Lgohtri-PBow TvAcc’'; 
end % End of Configuration I sub-Loop 


fe 


elseif channels==38 


1f£ ‘ch== 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch== 
fscaler 
psdunits 
raounits 
chtitle 
elseit ch== 
fscaler 
psdunits 
raounits 
chtitle 
elseif eh-—4 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch== 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch== 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch== 
fscaler 
psdunits 
Tracunits 
ehtrtle 
elseif ch== 
fscaler 
psdunits 
raounits 
chewede 


fsl; 
pus: 
ney hd 
‘Wave Ht Bow’ ; 


Pee 
pul; 
ee 
CSOnis Soule = 


1S DR 
pud; 
eZ 
TREES. Roi; 
fS2: 
pud; 
ne BA 
‘'TACS Pitch’; 
£S3- 

pus; 

naa pe ee 


ten 


£ atl 
aaa S 


ROL RES 


Law Re; 


CG ZACGC’ : 
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Pitch Rt’; 


elseif ch==9 
fscaler 
psdunits 
raounits 
ent me ke 
elseif ch==10 
fscaler 
psdunits 
raounits 
chtitle 
Slser e=ch=—i1 
fscaler 
psdunits 
raounits 
chtitle 
elseif ecn==12 
fsscaler 
psdunits 
raounits 
Cchertele 
elseif ch==13 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==14 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==15 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==16 
fscaler 
psdunits 
raounits 
Gheatle 
elseif ch==17 
fscaler 
psdunits 
raounits 
chtitle 


{s4; 
pu4; 
mua 


“TASS eeG ev ACe 


ie 
pud4; 
roe 
‘TACS 


£s4-; 
pu4; 
mu 


“THES (Ol ZACe 


=—Es4 ; 
pu4; 
Tale, 
PACS 


ES7- 
puz; 
ialize: 
‘Cher (Roi; 
fs 
Paz 
Buz 
‘Later akches 
ES3: 
joes 
sal S % 
eee 2 


fs3; 
pus 
rae, 
* Totes 


fs4; 
pu4; 
1g eee 
“Talkies 


7 


CG xAee 


Chae eet, 


Rei Rts 


CG ZAcc’ 


f 


f 


i 


Uj 


i 


PaubehwRiee 


/ 


elseif ch==18 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==19 
fscaler 
psdunits 
raounits 
chtitle 
elsentf-ch——70 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==21 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==22 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==23 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==24 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==25 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==26 
fscaler 
psdunits 
raounits 
chtitle 


fs4; 
pu4 ; 
Taelele 
ite 


fs4; 
pu4; 
teuliae 
eet r 


£sl1; 
pul; 
ig) abe 
PAGS 


i les 
iu 1? 
ao © el eae 
‘TACS 


Pods, 
pul 
Teles: 
STACS 


fs 1; 
pul; 
mu 
‘Load 


Esl. 
jene 
5g Oe 
‘Load 


fs4; 
pu4; 
ruse; 


‘BoomTip-HorizAcc’ ; 


fs4; 
pu4; 
rua; 


‘BoomTip-Vert Acc’; 


GG 


YACC =, 


CG XACC = 


LhtRelX’; 


LhtRelY’ ; 


LhtRelZ’; 


Mo Orthognl’; 


Mo 
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Parallel’; 


f 


f 


, 


f 


elseif ch==27 
fscaler 
psdunits 
raounits 
chtitle 
ersetmeen=—26 
fscaler 
psdunits 
raounits 
ehelele 
elserrech——29 
fscaler 
psdunits 
raounits 
EhEaele 
elseif ch==30 
fscaler 
psdunits 
raounits 
Ghee 
elsei@ en—=35 1 
fscaler 
psdunits 
raounits 
Chizisele 
elseif ch==32 
fscaler 
psdunits 
raounits 
CHETE le 
elseif "eh——33 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==34 
fscaler 
psdunits 
raounits 
chtitle 
elseif 2echn==35 
fscaler 
psdunits 
raounits 
Git tieve 


fs4; 
pu4; 
mele 


“LOW@iet— Seman VACC’ ; 


jee 
pu4; 
Wael 


’Lohtr4-Bow VAcc’; 


Es 
pu4; 
ran uill 


‘Lohtr3-Bow VAcc’; 


fs4; 
pu4; 
nails 


NOMNGrEZ -Bow VAGe = 


fs4; 
pu4d; 
Tae 


‘Lohtrl-Bow VAcc’ ; 


ESA 
pu4; 
nab UAL 


c 


fd 


a 


a 


‘Lohtr2-SMdShpVAc ’ 


fs4; 
pu4; 
uae 


‘’Lohtr2-PMdaShpVAc’ ; 


£s4-; 
pu4; 
bag ae 


‘Lohtr2=CELMdaSpvVAc’ ; 


fs4; 
pu4; 
mua 


‘EoOnesl-CLMaGSpVAc’ ; 


8] 


f 


a 


/ 


elseui ch==36 
fscaler 
psdunits 
raounits 
Ghtiete 
elseif ch==37 
fscaler 
psdunits 
raounits 
chtitle 
else 
fscaler 
psdunits 
raounits 
chtitle 
end % End of 


else 


[hea ele = — 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==2 
fscaler 
psdunits 
raounits 
chtitle 


elsert ch== 


fscaler 
psdunits 
raounits 
chtitle 


elseif ch== 


fscaler 
psdunits 
raounits 
chtitle 


aleeit -ch== 


fscaler 
psdunits 
raounits 
chtitle 


It 


fs4; 
pu4 ; 
rué4 > 


‘Lgohtr4-StrnLngAc’ 


fs4; 
pu4; 
Tamla, 


Wontra—-PStrni hee. 


rsa; 
pu4; 
has 


‘Lgohtr1-PBowTvAcc ’ 
SonEicoeen IIT sub- Bess 


fsi;> 
Due: 
vg 8 ed Dae 
‘Wave Ht Bow’ 


fsl; 
Du 
nb Alle 
“Sonix Sone’ 


fs2; 
pud ; 
wy) We 
TACS Rods” 3 


£s2; 
DU2; 
ngs Wag 
“TACS Pitch’ 3 
Ee5- 

bus; 

Pils: 


“TRCs Rolly Ro 
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elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
Glitedeewre 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 


fs3; 
pus; 
adie 
“TAGS 


£S3: 
Dust, 
1g ble) & 
TAGS 


fs4; 
pu4; 
1g) Be 
‘'TACS 


fs4; 
pu4; 
1g Us 
‘TACS 


10 


fs4; 
pu4; 
miele. 
‘"TACS 


ital 


Te Ss 
pu4; 
mae 
TAGS 


d2 


fs4; 
pu4; 
mua 
‘TACS 


13 


ES 
DuZ, 
rae 


Pateh: Rt’ - 


Yaw Rt’; 


CG ZAGCE’s 


CC YACc = 


CG XAcc’ ; 


CP ZACe- 


CE OVACG. 


“DD 5 65meno il: 


14 


£SZ: 
piu; 
aa baer 


‘DD. 7636 Pac. 
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eliseiifech==15 


fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
Chtcitele 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 
elseif ch==19 
fscaler 
psdunits 
raounits 
entatie 
elseif ch= 
fscaler 
psdunits 
raounits 
chtwtle 
elseif ch= 
Efscaler 
psdunits 
raounits 
chtitle 
elseif ch= 
fscaler 
psdunits 
raounits 
Cherelle 
elseif ch= 
fscaler 
psdunits 
raounits 
chtitle 


i 


fs4-; 
pu4; 
eA 
Die 55 CG ZACe = 


fs4; 

pu4; 

1G 

mph eI05 CG. LAece. 
1% 

isd: 

pu4; 

ee 

“OD 963 6 xAGcc. 
18 

rope ig 

pul, 

1 6 


“TACSI DDI6s= hela. 


fsl; 
pul; 
veld 


‘TACS1 DD963 RelY’; 


20 
tol; 
pul; 
nen ola le 


‘TAGS!’ DD963S RellZ ; 


oral 
fsl1; 
pul; 
ng, 


‘Lead Mo Orthogn! ; 


ORS, 
fsl1; 
bul; 
neg) gl 


‘Load Mo Parallel’; 


23 
fs4; 
pu4; 
Bud: 


‘BoomTip-HorizAcc’ ; 


84 


s 


f 


é 


f 


f 


f 


oe 


else 


fscaler = fs4; 
psdunits = pu4; 
raounits = ru4; 
chtitle =  Boomlup=vert Acc’ 


end % End of Cont auction ie Si weep 

end % End of the Overall Configuration Loop 

Base-lining the Raw, Time-based Data--------------------- 
x = data(:,ch); 


oe = xX-mean(xX); 


FFT Computation------------------~------------------------ 
4 = f£ft(x,n); 


Model-Scale PSD Computation------------------------------ 
PSD = ee cong (Y) in. 


Moving-window Averaging of Model-Scale PSDs for 
smoothing------------------------------------------------ 


halfptavg = ptavg/2; 
sum = Q; 
One ie = it DieawG 
sum = sum + PSD(k,1); 
if k > halfptavg 
psd(k-halfptavg,ch) = sum/k; 
else 
end 
end 
pts = window + halfptavg; 
ptavgplusone = ptavg + I; 
y = 1; 
for k = ptavgplusone:pts 
sum = sum + PSD(ptavg+j,1) - PSD(j,1); 
psd(k-halfptavg,ch) = sum/ptavg; 
J = lee ol 
end 
meansgr = mean(x.%*2) ; 
psdarea = (trap (@momeda,osat- -emy). 
psdscaler = meansqr/psdarea; 
psd(:,ch) = psda({(>,ch) *psdscaler- 
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fold 


Computation of Model-Scale RAOs-----------~--------------- 
< = pscle.2ch) ./psd (a 
mrao(:,ch) Soa.) > 


iH 


fo) 


Computation of Full-Scale RAOs--------------------------- 
frao(:,ch) = mrao(:,ch)*fscaler: 


fe\e) 


Plot Model-Scale PSDs and Save Figure in jpg Format------ 
elt 
Diet {momega Cl: jose) er) ) 


Gils ee 
sch =seMidimescer (Ch) = 
title([’Model-Scale PSDs for DTMB Run ‘srun,’ / Channel 


‘seh 2) Cheaele | } 
xlabel(’Omega (rad/s) ’) 
ylabel([’PSD ‘psdunits] ) 


saveas(gcf, [swrite,folder,’\’,filename,’C’,sch,'’p’],'’jpg 
‘); & D:\###\REFECHED. JPG 


oe 


Plot Model-Scale RAOs and Save Figure in jpg Format------ 
elie 


plot (momega(l1,:),mrao(:,ch)) 

grid 

title([’Model-Scale RAOs for DTMB Run ‘Srun,’ / Channel 
(sen, 3: “ehtrele)} 


xlabel(’Omega (rad/s) ') 
ylabel([’‘RAO ‘raounits] ) 


saveas(gcf, [swrite,folder,’\’,filename, ’C’,sch,’m’],’Jpg 
‘); ©& D:\###\REF#FCHEM. JDG 


fel 


Plot Full-Scale RAOs and Save Figure in jpg Format------- 
clf 


DloOt (fomedgai(i,:) .frac{:, ch) + 

grad 

title([’Full-Scale RAOs for DTMB Run ‘’srun,’ / Channel 
are) ol a ee © os al 


xlabel(’Omega (rad/s) ’) 
ylabel([{’RAO ‘raounits] ) 


saveas(gcf, [swrite, folder, ’\’‘,filename, ‘C’,sch,'f£’'],'jJpg 
‘); %& D:\###\REFFCHEE. JDO 


end % End of all Channel Looping 
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% Save Omega Matrices as Last Column of Composite PSD and 
RAO Matrices------------------------------------+-~-------- 


ela! —aelglc 1b 5 

psd (=: 2en) = momega(l1,:)'’; 
Mascl@aG Cl)) = momega(l,:)’; 
ERAOKL en) = fomega(l,:)'; 


Save Composite PSD and RAO Matrices---------------------- 
D: \###\R###psd. txt 
D:\###\R###mrao. txt 
D: \###\R###£rao.txt 


0d dP oP 


oe 


save ([{swrite,folder,’\’,filename, '’psd.txt’],’psd’,’-ascii 
mals vas 

save ([swrite, folder,'’\’,filename, ‘mrao.txt’],’mrao’,’-ascii 
in Stabs ) = 

save([swrite,folder,’\',filename,’frao.txt’],’frao’,’-ascii 
‘,’-tabs'); 


th TGTONOL Ge TRUM! VEN gt SYS per 


° 
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ie) 


ol ke 


APPENDIX H. PTAVG COMPARISON PROGRAM 


% Program title: ptavgcomparison.m------------------------- 
clear 


Load Spreadsheet Results for run 67---------------------- 


load 


(e) 


load 
load 
load 
load 


load 
load 
load 
load 


load 
load 
load 
load 


%* Load RAO Program Results for ptavg = 


Iai Owes excl ° 

50, 76 and 100-- 
ROUGHWDSaZ6 .c xe - 

ROO /psd50. EXE; 

RO67psda76.txt; 

R06 pseal00 Sexe - 


ROG /mraoco. exe; 
RUG/MrEaAOsU exe: 
ROG6G/Mrac/ 6 —eExE- 
R06 /meaol00. EXE; 


ROG 7ECaAOZ6 Cx; 
RUG? Braos 0] Gice; 
RUGy £rao76. exe. 
RUG? Eraol00- EXE; 


% Plot Model-Scale PSD Comparison-------------------------- 

figure (1) 

plot (RUG 7pSaz26(>, 4/7) RUG /psd261-7 3 eae. 
RO67psd50(:,47),RO67psd50(:,3),... 
RO67psd76(:,47),RO067psd76(:,3),... 
Rec psalo0 (: 47 ROG psalte(:.3) ou... 

mamMow |. 2 )e, uno? (52,4 )5} 

title(’RO67CH3 Model-Scale TACS Roll PSDs 

Comparison’ ) 

xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (deg) ’) 

grid 


iu ptavg iT 
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% Plot Model-Scale RAO Comparison-------------------------- 

figure (2) 

plOGE(RO6 /mraoZz6(: ,47) Rv e mr a0Z26 (=, 3) 
RO6/7mrao50(:, 47 )@RU@e7mracs0 (: 3). 42. 
ROG /mraco/76(:,47)RUG7mrao/6(:, : 
RO6G/mraol00(: 47) Res /mrachoee 3) >). 

Puno / (2,2), PUnG7 (2 See) 

title(’RO67CH3 Model-Scale TACS Roll RAOs "ptavg" 
Comparison’ ) 

xlabel(’Omega (rad/s) ’) 

ylabel(’RAO (deg/in) ’) 

Grid 


% Plot Full-Scale RAO Comparison--------------------------- 
figure (3) 
DlOtU( R06 frae7G 440) Reo PraoZz6t. 3) pe 
ROG JEraOSU sea 7) ORUG (i rao50(: 3) 3a 
ROG? Eaes 61047) ~RUC Erao 1/6) eee 
RO67Eracl0C (47 PROG 7 Eraol00 (3) 
title(’RO67CH3 Full-Scale TACS Roll RAOs "ptavg" 
Comparison’ ) 
xlabel(’Omega (rad/s) ’) 
ylabel(’RAO (deg/in) ’) 
Giac 


% End of ptavgcomparison program----<-<------<<--<-<-=<-=sse= 
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APPENDIX I. PSD AND RAO COMPARISONS PROGRAM 


Program Title: psdandraocomparisons .m-------------------- 


Computational and Graphical Comparison between the 
Various Met Heys —— = a ee 


clear 
cle 


format short 


° 
© 


load ROQ67; 

beac HWinNG 7 4 EX ; 
load R067psd. txt; 
load RO67mrao.txt; 
load RO67frao.txt; 


° 
© 


Load Applicable Matrices--------------------------------- 
Raw, Time-based Data 

Spreadsheet PSD and RAO Results 

RAO Program Model-Scale PSD Results 
RAO Program Model-Scale RAO Results 
RAO Program Full-Scale RAO Results 


of bo cP oO 


fo) 


Base-— |“imeomine Raw emime—-based—-—— ] — = a ee eee 
=. Gata) % Wave Ht Data 
= X - mean(X); 


= eClel ee.) % TACS Roll Data 
= Z - mean(z); 


Computation of Wave Ht PSDs and Associated Frequencies--- 


[hwelch, fwelch] = -pwelcin@e, [leh ,s0g6, 3242): 

[hperio, fperio] = periodogram(x, [],8096,32.2); 

tamales, Fnuita| = oueincs Ne lS Gas 2", As 

meansqrx = mean(x.*2); % Wave Ht Data 
welchareah = trap (wwelch,hwelch) ; % Wave Ht 
perioareah = trap(wperio,hperio) ; 

multiareah = Gta (wml et nim ip ye 

hwelchscaled = hwelch*meansgqrx/welchareah; 
hperioscaled = hperio*meansgqrx/perioareah; 
hmultiscaled = hmulti*meansgqrx/multiareah; 


9] 


6 Computation of TACS Roll PSDs and Associated Frequencies- 


[rwelch, fwelch] = pwelch(z,[],[],8096,32.2); 

[rperio, fperio] = periodogram(z,[],8096,32.2); 

Baiulbicay aubiucam|, = pinta | 1) 6096,32.2): 

meansqrz = mean(z.%2); % TACS Roll Data 
welcharear = trap(wwelch,rwelch) ; *® TACS Roll 
periloarear = trap(wperio,rperio); 

multiarear = Crap (wm ti, rmulta)- 

rwelchscaled = rwelch*meansgqrz/welcharear; 
rperioscaled = rperio*meansgrz/perioarear; 
rmultiscaled = rmulti*meansgqrz/multiarear; 


6 Convert Frequency to Omega------------------------------- 


wwelch =e Oa hwo en. 
wperl1o = 2*pi*fperio; 
wmulti =e een vague cae. 


% Compute Areas under the Wave Ht PSD Curves---------——— 2 


welchareah2 = trap(wwelch,hwelchscaled) ; 
perioareah2 = trap(wperio, hperioscaled) ; 
multiareah2 = trap(wmulti,hmultiscaled) ; 
spreadareah =. trap (runes (: 72), 2uUn6 7 ees) 
raoprogareah = trap R067psa - 64 RUG eset - a 


% Compare Areas under the Wave Ht PSD Curves against 
meansqrs-—------ oe er rrr rr ee 
wavehtpsdareacomparison = [meansgqrx welchareah2 
perioareah2 multiareah2 spreadareah raoprogareah] 


%¢ Compute Areas under the TACS Roll PSD Curves---------==== 


welcharear2 = trap(wwelch, rwelchscaled) ; 
perlioarear2 = trap(wperio, rperioscaled) ; 
multiarear2 = trap(wmulti,rmultiscaled) ; 
spreadarear =“trapirun67{>,2), seme) (274) )- 
raoprogarear = traptR06/psd(:, 47), R067 psa :,.3)8 


% Compare Areas under the TACS Roll PSD Curves against 
meansqrs------------------------------------------------- 


rollpsdareacomparison = (meansgqrz welcharear2 
perioarear2 multiarear2 spreadarear raoprogarear] 


22 


% Graphical Model-Scale Wave Ht PSDs Comparison of All 5 
Methods------------------------------------- -- -- -- 


figure (1) 

plot (wwelch, hwelchscaled) % Welch Method 

nota. On 

plot (wperio, hperioscaled) % Periodogram Method 
plot (wmulti, hmultiscaled) % Multi-taper Method 
muerte / (2,2), cuUn6) (2,3) % Spreadsheet Results 
mlGe (ROO 7pSal 47) RUGIosa(=. )))- % Rao Program Results 


title(’RO67CHO1 Model-Scale Wave Ht PSDs Comparison’ ) 
xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (dimesionless) ’) 

eecrs( (0 6 0 Oda» 

Chenke| 


% Graphical Model-Scale Wave Ht PSDs Comparison------------ 
Spreadsheet vs. RAO Program only 


oe 


figure (2) 

Pao t (eG) (ae. 2), HUM 6 7 (2. 3).) % Spreadsheet Results 
Bone) On: 

meOe(RUG/psa(., 48) ,;ROUGJpsa(-<-. ))) % Rao Program Results 


title(’RO67CHO1 Model-Scale Wave Ht PSDs Comparison’ ) 
xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (dimesionless) ’) 

erates (1096 0 20065] } 

eiaicl 


% Graphical Model-Scale TACS Roll PSDs Comparison of All 5 
Miele la eChlee = = SSS SS = = oS SS oS 5 5 5 5 SS 5 SS SS SS eS 


mrgure (3) 

plot (wwelch, rwelchscaled) % Welch Method 

me la. on 

plot (wperio, rperioscaled) % Periodogram Method 
plot (wmulti, rmultiscaled) % Multi-taper Method 
Pik@ Gino? (2 2 eiaminG 7 (4) ) % Spreadsheet Results 
mor {( R06 /psa(7,4/), ROG7psat: 75.58 % Rao Program Results 


title(’RO67CH03 Model-Scale TACS Roll PSDs Comparison’ ) 
xlabel(’Omega (rad/s) '’) 

ylabel(’PSD (deg) ’) 

exsis [0-76 0°. 02)) 

Speake 


2, 


6 Graphical TACS Roll PSDs Comparison---~-~----------------- 


% Spreadsheet vs. RAO Program only 

figure (4) 

PEomieunG) (:,2) ,/ runes) % Spreadsheet Results 
held on 

PEOELRIUG7pSal: 47) 7 RUG esa, 3)) % Rao Program Results 


title(’RO67CH03 Model-Scale TACS Roll PSDs Comparison’ ) 
xlabel(’Omega (rad/s) ’) 

ylabel(’PSD (deg) ’) 

axes ( (UNGeG9 2002) 

Grad 


% Compute Model-Scale Roll RAOs from MATLAB Methods-------- 


R1 = rwelchscaled./hwelchscaled; 
welchmrao =" JS Cran Remi 
R2 = rperioscaled./hperioscaled; 
periomrao = Shit hee. 
R3 = rmultiscaled. /hmultiscaled; 
multimrao =" “Seite 1 


% Graphical Model-Scale TACS Roll RAOs Comparison of All 5 
Methods--—--—-=---- = — == = = = = = 5 = 8 5 8 5 = = = = 


figure(5) 

plot (wwelch, welchmrao) % Welch Method 

hold on 

plot (wperio,periomrao) % Periodogram Method 
plot (wmulti,multimrao) % Multi-taper Method 
DE@elwunes (272) ne 7 Ot}. 5) % Spreadsheet Results 
pIoE (ROG 7mrac(:,47),R06/mrao(:,3)); % Rao Program Results 


title(’RO67CH03 Model-Scale TACS Roll RAOs Comparison’ ) 
xXxlabel(’‘Omega (rad/s)') 

ylabel(’RAO (deg/in) ’) 

axis([{0 6 0 10]) 

Gira 


% Graphical Model-Scale TACS Roll RAOs Comparison----------- 


% (Spreadsheet vs. RAO Program only) 

figure (6) 

plot (rune 7 (:.2) ,funoe? {:35)) % Spreadsheet Results 
neta. on 

plot (RO67mrao(:,47),RO6/mrao(:,3)) % Rao Program Results 
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title(’RO67CHO3 Model-Scale TACS Roll RAOs Comparison’ ) 
xlabel(’Omega (rad/s) ’) 

ylabel(’RAO (deg/in) ’) 

peas (10) 6 01252) ) 

Grid 


* sCale to Full-Scale RAOs and Omegas---------~------------ 


welchfrao = welchmrae*12/7/24.175:; 
periofrao = Der OMmsag~ 7424 do 
multifrao Seemed inisae ee 27 2h) oe 
wwelchfull = wwelen/saqrt (24.175) 
wperiofull = Whe de, = Cian 240 a 
worl ta ful) =" WMUl ti / Samet 24. le >). 


% Graphical Full-Scale TACS Roll RAOs Comparison----------- 


% (4 Methods Not including Spreadsheet) 
figure(7) 
plot (wwelchfull,welchfrao) % Welch Method 


hele! on 

plot (wperiofull,periofrao) 

plot (wmultifull,multifrao) 

MPuOEmrrt R06 Jirrao(:, 47), R0G/frao(:, 3) } 


oo 


Periodogram Method 
Multi-taper Method 
RAO Program Results 


oe 


feed 


title(’RO67CHO3 Full-Scale TACS Roll RAOs Comparison’ ) 
xlabel(’Omega (rad/s) ’) 

ylabel(’RAO (deg/in) ’) 

peeren | Owl. 0. 5 |) 

Ghagikel 


6 Graphical Full-Scale TACS Roll RAOs Presentation--------- 
RAO Program only 

figure (8) 

MbOt (UG) Erao(:), 47 ye RUCTErao( 5 )) % Rao Program Results 


oP 


title(’RO67CHO0O3 Full-Scale TACS Roll RAOs Presentation’ ) 
xlabel(’Omega (rad/s) ’) 

ylabel(’RAO (deg/in) ’) 

eeaicreiO 1-4 0 1.47} 

opgalye| 


% End of PSD and RAO Comparison Program-------------------- 
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